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La presente invention a pour objet une methode pour deformer graduellement 
une repartition initiale d'objets de nature geologique, formee par simulation d'un 
mod^e stochastique de type objet, d'aprfes des mesures ou observations, pour 
r adapter au mieux k des contraintes physiques impos^es, de type hydrodynamique 
5 par exemple. 

La methode selon T invention trouve des applications dans la modelisation de 
zones souterraines o\X il s'agit de gendrer des representations montrant comment est 
distribufe une certaine grandeur physique dans une zone du sous-sol (la perm6abilite 
par exemple), qui soient compatibles au mieux avec des donnees observees ou 
10 mesurees : donnees geologiques, enregistrements sismiques, mesures obtenues dans 
■ des puits notamment mesures des variations au cours du temps de la pression et du 
debit de fluides issus d'un gisement, etc. 

ETAT DE LA TEGHNIQUE 

Par le brevet FR 2 780 798 du demandeur, on connait une methode pour 
15 deformer graduellement un modele stochastique (de type gaussien ou appareiite) d'un 
milieu heterogene tel qu^une zone souterraine, contraint par un ensemble de 
parametres relatifs a la structure du milieu. Cette methode comporte le tirage d*un - 
nombre p de realisations (ou representations) independantes du modele ou d'une 
partie au moins du modele choisi du milieu, parmi I'ensemble de toutes les 
20 realisations possibles et une ou plusieurs etapes iteratives de deformation graduelle 
du module en effectuant une ou plusieurs combinaisons lineaires successives de p 
realisations initiales independantes entre elles puis des realisations composites 
successivement obtenues 6ventuellement avec de nouveaux tirages, etc., les 
coefficients de cette combinaison etant tels que la sonmie de leurs carres est egale h. 
25 1. 

Par le brevet FR 2 795 841 du demandeur, on connait une autre m6thode pour 
deformer graduellement les representations ou realisations, generees par simulation 
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sequentielle, d'un module stochastique non ndcessairement gaussien d'une grandeur 
physique z dans un milieu heterogene maill6, afin de les ajuster k un ensemble de 
donnees relatives h la structure ou T^tat du milieu qui sont collectdes par des mesures 
et observations prealables. La m6thode comporte essentiellement Tapplication d'un 

5 algorithme de deformation graduelle d'un modele stochastique a un vecteur gaussien 
a N variables mutuellement independantes, qui est relie a un vecteur uniforme a N 
variables uniformes mutuellement independantes par la fonction de repartition 
gaussienne de fagon a defmir des realisations du vecteur uniforme, et I'utilisation de 
ces realisations pour g6nSrer des representations de cette grandeur physique z, que 

10 Ton cale par rapport aux donnees. 

Les m^thodes pr6c6dentes sont applicables aux modeles mailles (modeles de 
type pixel) qui conviennent pour modeliser des champs de grandeurs continues et, de 
ce fait, sont mal adaptes h la modelisation de zones traversees par des rdseaux de 
fractures ou des systemes de chenaux, par exemple. 

15 Les modeles bases sur des objets sont des arrangements dans Tespace, d'une 

population d'objets definis g6ometriquement. A la base, un modele de type objet est 
un module bool6en que Ton peut d6finir conmie une reunion d' objets identiques par 
nature, distribuds au hasard dans Tespace. Les modeles bool6ens (de type objet) 
presentent un gr^d interSt pour la description g€om6trique des milieux heterogenes, 

20 tels que les systemes de depot m^andriformes, les reseaux de fractures, les milieux 
poreux ^ Techelle granulometrique, les milieux vacuolaires, etc. Les objets 
geologiques sont definis par leurs formes et leurs tailles. Leurs emplacements dans le 
champ sont definis en tenant compte de leurs interactions : attirance-repulsion, 
tendance au regroupement (clustering), etc. 

25 Contrairement aux modeles de type pixel, les modeles bases sur des objets 

peuvent foumir par exemple des representations geologiques realistes d'un reservoir 
souterrain a un stade precoce ou les donnees obtenues par des mesures in situ sont 
encore rares. 
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L*6tat de Tart dans le domaine des modeles de type objet, est decrit 
notamment dans les publications suivantes : 

- Matheron, G., 1967, "Element pour une theorie des milieux poreux", Masson, 
Paris; 

5 - Matheron, G., 1975, "Random sets and integral geometry", Wiley, New York; 

- Serra, J., 1982, "Image analysis and mathematical geology", Vol.1, Academic 
Press, London; 

- Stoyan, D.S., et al,, 1995, "Stochastic geometry and its applications", 2*^** Edition, 
Wiley, Chichester; 

10 - Lantuejoul, C, 1997, Iterative algorithms for conditional simulations, in Baafi et 
others, eds.; Geostatistics WoUongong 96, VoL 1, Kluwer Acad. Publ., 
Dordrecht, The Netherlands, p.27-40. 

La position des objets dans un modele de type objet est distribute suivant un 
processus ponctuel de Poisson bien connu des specialistes. La forme et la taille des 
15 objets sont independantes de leurs positions. Ce mod&le peut etre generalise par une 
combinaison d' objets de nature diff6rente ou/et en utilisant un processus ponctuel de 
Poisson de densite non-stationnaire. i 

Bien que les modeles booleens aient 6i6 largement etudies dans la litterature, 
on ne dispose pas de m6thode coherente et efficace pour contraindre ces modeles aux 

20 donnees physiques, notamment aux donnees hydrodynamiques, ce qui est pourtant un 
enjeu majeur pour leur application en ingenierie de reservoir. Les methodes de 
deformation graduelle des modeles stochastiques de reservoir de type pixel telles que 
d6crites par exemple dans les deux brevets precites, ne sont pas directement 
utilisables pour les moddles booleens. Contraindre les modeles booleens aux donnees 

25 hydrodynamiques par exemple requiert le developpement des algorithmes coherents 
pour la deformation et le deplacement des objets. 
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LA METHODE SELON L'lNVENTION 

La methode selon Tinvention permet de generaliser la technique de 
deformation graduelle ddcrite dans les deux brevets precites, aux modeles booleens 
stationnaires ou non-stationnaires et ceci avec ou sans contrainte geometrique aux 
puits. La methode s'avere particulierement utile notamment pour les ingenieurs de 
reservoir soucieux d'ajuster de fagon coherente et efficace les modeles de reservoir 
de type objet. 

La methode selon Tinvention permet de deforaier graduellement une 
realisation initiale ddfinissant la repartition d'un ensemble d'objets dans une zone 
d*un milieu hdterog&ne tel qu'une structure geologique, formee par simulation d'un 
modfele stochastique de type objet, les objets etant repartis dans la zone selon un 
processus ponctuel de Poisson sous forme de points figuratifs avec une densite de 
points X{x) qui varie en fonction de leur position (x) dans la zone. EUe comporte 
essentiellement les deux etapes suivantes : 

- on genere une realisation d'un vecteur aleatoire uniforme selon lequel la position 
de chaque objet est definie tout en respectant la densite A(x) ; et 

- on modifie graduellement le vecteur aleatoire uniforme selon une procedure de 
defomiation graduelle, pour obtenir la migration graduelle de chaque objet et par 
consequent le changement graduel de la repartition des objets dans la zone, 
jusqu'^ obtenir une realisation finale ajust^e au mieux a des paramStres relatifs ^ 
la structure du milieu tels que des parametres hydrodynamiques. 

n est possible de limiter la migration des points figuratifs d'objets dans un 
sous-domaine de la zone (un puits traversant la zone par exemple) en imposaht une 
densite de points nuUe dans la partie complementaire du sous-domaine. 

Suivant un mode de mise en ceuvre, on passe graduellement d'une realisation 
contenant un premier ensemble de TV, points ^ une realisation contenant un 
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deuxieme ensemble de N2 points en construisant une chaine N{t) de nombres de 
Poisson entre les deux nombres et en utilisant la procedure de d6foraiation 
graduelle. 

H est possible aussi de modifier graduellement la taille, la forme et 
Torientation d'un objet au cours de sa migration en utilisant la procedure de 
deformation graduelle. 

En cas o\l elle est entach^e d'incertitude, il est possible aussi d'ajuster 
graduellement la densit6 de points A(jc) en utilisant la procedure de deformation 
graduelle. 

PRESENTATION SOMMAIRE DES FIGURES 

D'autres caracteristiques et avantages de la m6thode selon Tinvention, 
apparaitront h la lecture de la description d6taillee ci-apres, en se refdrant aux dessins 
annexes oCi : 

- la Fig.lA, IB montrent respectivement une fonction de densite d'un processus 
ponctuel de Poisson non-stationnaire et une realisation du processus ponctuel de 
Poisson non-stationnaire g^n^r^e par la m6thode sdquentielle ; 

- les Fig.2A a 2D montrent differents exemples de' trajectoires de migration 
graduelle entre deux points d*un processus ponctuel de Poisson non-stationnaire ; 

• les Fig.3A a 3H montrent diff^rentes etapes successives d'une realisation du 
processus ponctuel de Poisson non-stationnaire lors de la migration graduelle des 
- points ; 

- la Fig.4 montre les domaines possibles de migration d'un disque dans le cas de 
trois points conditionnant ; 
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- la Fig.5 montre a titie d*exemple, une chaine complete de realisations successives 
d'une simulation bool6enne d'objets elliptiques ou le parametre de d6foraiation t 
varie entre - tt et ;r par pas At =0.1 tt ; et 

- la Fig.6 montre a titre d'exemple, une chaine incomplete de realisations 
5 successives d'une simulation booleenne d'objets elliptiques oH le parametre de 

deformation t varie entre 0 et 0.2 ;r par pas =0.01 n . 

DESCRIPTION DETAILLEE 
Generalites 

Les objets geologiques auxquels la methode s' applique sont par exemple des 
10 fractures a plus ou moins grande echelle, ^ I'interieur d'une zone reservoir, ou bien 
encore des chenaux. La methode peut s'appliquer egalement a des structures 
granulaires ou vacuolaires de taille beaucoup plus reduite. Tous ces objets sont 
difficiles k mod61iser par des modeles de type pixel. 

Les operations de migration progressive qui vont etre decrites, s'appliquent a 
15 un modele initial oil les positions des objets sont repr6sent6s par des configurations 
de points (dites processus ponctuels). La repartition de ces points varie en densite en 
fonction de leurs positions dans la zone modeiis^e. Cette repartition est basee sur 
differentes donn6es connues par mesures ou observations : mesures gdomecaniques 
obtenues dans des puits par exemple, donnees sismiques issues d' operations 
20 sismiques prealables. 

Partant de ce modele initial, la methode va permettre de deformer 
graduellement la repartition initiale suivant une serie de regies, de maniere que dans 
la repartition finale, le modele soit optimise pour mieux correspondre a des 
contraintes que Ton impose sur un ou plusieurs parametres physiques, telle que par 
25 exemple une distribution de valeurs de permeabilite. Les regies de migration sont 



1er depot 



telles que Ton peut deplacer globalement un grand nombre de points diff6rents du 
moddle k partir d'un nombre reduit de parametres de controle. 

En general, on definit une fonction objectif qui mesure la difference entre les 
parametres physiques issues du milieu h6t6rog6ne reel et ceux simul^s sur une 
5 r6alisation du moddle stochastique. La valeur de la fonction objectif depend done des 
parametres de controle du modfele stochastique. On obtient les valeurs de ces 
parametres de contrdle par minimisation de la fonction objectif. 

Rappel sur le processus ponctuel de Poisson 

Le processus ponctuel de Poisson est un ensemble aleatoire denombrable de 
10 points r^partis dans tout espace SH" . Cet ensemble de points a les caract^ristiques 
suivantes: 

Soit D un domaine de SH" . Si le volume de D , note |d| , est fini, alors le 
nombre de points tombes dans D , note N(D) , suit une loi de Poisson de paranietre 
A|D|.Soit 

15 p[iV(D) = n]=e-^>^^^J^ Vn^O (1) 

n! 

ot A est appelee la densite du processus ponctuel; elle mesure le nombre moyen de 
points tombes dans un domaine de volume unitaire de Si'* , 

Soient Z)i,D2,...,£>^ des domaines de 91" disjoints deux a deux, alors les 
nombres de points tombes dans ces domaines iV(£>i),iV(D2 ),-.., ^(£^jk) sont des 
20 variables aleatoires mutuellement inddpendants. 

Conditionnellement a iV(D) = np, ces points sont independants et 
unifomiement repartis dans D . 
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Dans ce qui suit, nous nous interessons au processus ponctuel de Poisson dans 
le domaine borne D . 

Migration d'un processus ponctuel de Poisson stationnaire 

Considerons le probleme de migration d'une realisation d'un processus 
5 ponctuel de Poisson stationnaire dans D rectangulaire. Afin d'alleger la presentation, 
on admet que D est un hyper cube unitaire [0,1]" a n dimensions. Soientjc, et 
deux points ind6pendants unifonn€ment tires dans [0,1]" . Nous definissons une 
trajectoire entre et 



10 ou G est la fonction de repartition gaussienne centric et reduite. Selon Talgorithme 
de deformation graduelle decrit dans le brevet FR 2 780 798 precite, pour tout t , 
x{t) est un point uniforme dans [0,1]". Quand les deux points sont fixes, la 
trajectoire de la migration graduelle entre eux est entierement d6termin6e. Le 
changement d'emplacement d'un des deux points changera la trajectoire de migration. 

15 On peut montrer que la trajectoire definie par I'equation (2), est sym€trique par 
rapport au centre du domaine [0,1]" , et ceci quel que soit le nombre de dimensions 
n. Cela suggere que meme si les deux points sont isdes dans un coin/cote du 
domaine, la trajectoire de la migration graduelle entre eux peut toujours atteindre la 
partie opposee du domaine. 

20 Migration d'un processus ponctuel de Poisson non-stationnaire 

On considere maintenant le processus ponctuel de Poisson dans un domaine 
D de densite X{x) de fomie generale. Le nombre de points dans D est une variable 

aleatoire de Poisson de moyenne X{D) = ^^?i{x)dx. Ces points sont distribues 

independamment dans D suivant la fonction de densite de probability : 



jc(0 = G[G~H^i)cosr + G"*(jc2)sin/] 



(2) 



1 er depot 




9 



f{x)^k{x)lX{D) xeD (3) 

La simulation d'un processus ponctuel de Poisson de densite X(x) dans D 
peut se realiser en deux 6tapes: 

- generer un nombre n selon la loi de Poisson de moyenne egale a A(Z)) , puis 

- generer n points dans D independamment les uns des autres selon la meme 
densite de probabilite f{x) . 

On peut voir sur Texemple de la FigJA une fonction de densite, et sur celui 
de la Fig. IB, une realisation de processus ponctuel de Poisson de densite A(x) , 

Si on simule la loi f(x) par I'inversion de sa fonction de repartition, le point 
X correspond alors a un vecteur uniforme u . On peut done appliquer Talgorithme de 
migration graduelle au processus ponctuel de Poisson de density Mx) . Par 
construction, cette methode preserve la density et le nombre de points du processus 
initial. Les Fig,2A k 2D montrent quatre exemples de trajectoire de migration, 

•J. 

Prenons Texemple du processus ponctuel de Poisson non-stationnaire dans un 
domaine en deux dimensions D=[0,l]^ dont la densite est lin^airement croissante 
dans Taxe X et constante dans Taxe Y. Soit (x,y) le vecteur des coordonnees d'un 
point dans D de ce processus. Alors (x^y) admet la densite de probabilite 
bivariable : 

fix, y) = 2x, (x, y) e [0,1]\ (4) 

La simulation d'un point selon la loi (4) precedente est simple: 

- generer I'abscisse x selon la densite lineaire f(x) = 2x , puis 

- generer I'ordonnee y uniformdment entre 0 et L Soient 
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(5) 



y = v 

ou w et V sont deux nombres independants et uniformes entre 0 et 1. 

Nous pouvons ainsi appliquer ralgorithme de migration graduelle (2) au 
vecteur uniforme (m,v) afin d'etablir une trajectoire de migration du point (x,y) 
dans D. 

En pratique, la fonction de density A(jc) est souvent sous forme de grille. A 
titre d'exemple, on considere le cas d'une grille en 2 dimensions de M x N noeuds. 
Soient et yj les coordonnees du noeud (z\ j) . La loi marginale de x^ est alors : 



A partir de ces lois de probabilite monovariables, il est aise de generer les 
points du processus de Poisson de densite M^) sur une grille. 

Afin d'illustration, nous avons construit une fonction de densite en deux 
dimensions a partir d'une simulation gaussienne centree et reduite. Le variogramme 
est anisotrope et de type gaussien. Les directions principales d' anisotropic sont 
diagonales par rapport au systeme des coordonnees. Les facteurs d'dchelle dans les 
directions principales d' anisotropic sont de 0.3 et 0.1 respectivement. La taille du 
champ est de 1x1 et il est discrdtise en 1000x1000 pixels. Les nombres gaussiens 
sont transformes en nombres positifs selon T expression suivante : 



/U) = £/U.>>'y) (6) 



et la loi conditionnelle de yj sachant x^ vaut : 



f.Syj)^f(^nyjyn^i) (7) 



A(jc) = 400e^'<'^ 



(8) 



oil y(jc) est la simulation gaussienne. La Fig. lA montre la fonction de densite ainsi 
construite. La figure IB montre une realisation de processus ponctuel de Poisson 
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admettant la fonction de densite de la Fig lA. 2000 points sont generes selon la 
mdthode s6quentielle. 

Les Fig.3A k 3H illustrent revolution d'une realisation du processus de 
Poisson lors de la migration graduelle des points. On peut constater qu'au cours de la 
migration, la densite de points est respect6e. 

Migration d'objets de simulations booleennes conditionnelles 

La methode ci-dessus peut etre immediatement appliquee a la migration d'un 
domaine S de forme quelconque dans D . En fait, la migration d'un point dans S peut 
etre realis6e dans D en utilisant la fonction de density de probabilite tronquee : 

/5U) = A(;c)1^5/A(5), xeD, 5cD. (9) 

Par I'algorithme d*it6ration markovienne decrit dans la reference k 
Lantuejoul, 1997, on peut simuler un modele booleen dans un domaine D sachant 
que deux sous-ensembles Q et Q de D appartiennent respectivement a T union des 
objets et a son complement. Ensuite, les algorithmes de migration dans un domaine 
quelconque peuvent 6tre utilises pour la deformation graduelle des simulations 
booleennes conditionnees par des donnees g^ometriques des puits. En effet, a partir 
d'une simulation booleenne conditionnelle et sans compromettre le conditionnement 
par Q et les objets doivent se deplacer uniquement dans leurs domaines 
respectifs d€finis selon les formes des objets et la configuration de C, et Cq - 

Considerons un objet A d'une realisation bool6enne conditionnelle, qui 
inclut un sous-ensemble de Q et exclut Q - Si lors de la migration de A , il doit 
toujours inclut mais exclut Q , le domaine autorise de migration de Tobjet A 
est: 

D^={x:ConA,=0;Q^cAj (10) 
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La Fig. 4 montre un cas avec trois points conditionnant et les huit domaines 
possibles de migration d*un disque. Si par exemple le disque doit toujours couvrir les 
points (a) et (b) mais eviter le point (c), alors son centre peut se d€placer seulement 
dans le domaine 3. 

Migration avec apparition et disparition de points 

Le nombre de points dans D d'un processus ponctuel suit une loi de Poisson 
de parametre A(D) . H est done n6cessaire de faire varier le nombre de points dans D 
au cours de leur migration. Dans cette section, nous presentons d'abord deux 
methodes pour constmire des chaines de nombres poissonniens, et ensuite nous 
verrons comment migrer entre deux ensembles de points dont les cardinaux ne sont 
pas identiques. 

Nous cherchons k constmire une chaine de nombres poissoniens entre deux 
nombres poissoniens N^ et TV^* generes independamment par Tinversion de la 
fonction de repartition. Soient et U2 deux nombres uniformes (entre 0 et 1) 
independants k partir desquels sont obtenus les nombres et 



ou F"* designe la fonction de repartition inverse de la loi de Poisson. Selon 
Talgorithme de deformation graduelle, nous pouvons construire une chaine de 
nombres uniformes entre Uj^ et U2 par 



En inversant la fonction de repartition de la loi de Poisson, on obtient une 
chaine de nombres poissonniens 



(11) 



U{t) = G[G-' (U, ) cosf + G'' ) sin t] . 



(12) 



Nit) = F-'lU(t)]. (13) 
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Le calcul de la fonction de repartition inverse peut s'effectuer par dichotomie. 
N^anmoins, si le paramfetre de la loi de Poisson est trop grand, cette methode reste 
quelque peu couteux. 

Pour 6viter le calcul de la fonction de repartition inverse de la loi de Poisson, 
5 on peut envisager une autre fa5on de g^n^rer les nombres poissoniens. On sait que le 
nombre de sauts, d'un processus de Poisson de param^tre unite, dans un segment de 
longueur A, suit pr^cisement une loi de Poisson de paramStre A. Simuler un 
processus de Poisson est chose simple. En effet, les intervalles entre deux sauts 
consecutifs du processus sont independants et suivent la loi exponentielle de 
10 parametre 1. En deformant graduellement les segments exponentiels, on obtient une 
chaine de nombres poissoniens. H est aise de construire une chaine de nombres 
exponentiels, car le calcul de la fonction de repartition inverse de la loi exponentielle 
est simple. On rappelle ici que la fonction de repartition de la loi exponentielle de 
paramdtre 1 s'dcrit 

15 F(s) = l-e-\ s>0 (14) 

et sa fonction inverse s'ecrit : 

F^' (u) = -ln(l - u), 0 < «< 1 . (15) 

Toujours est-il que, conrnie le nombre de sauts dans le segment de longueur 
A est en moyenne egal a A , le nombre de chaines de nombres exponentiels est en 
20 moyenne proche de A . 

Considerons maintenant la migration d'un ensemble de points (ensemble 
1) k un ensemble de N2 points (ensemble 2). Conmie ^N^^ la migration de 
Tensemble 1 k I'ensemble 2 implique n^cessairement Tapparition ou la disparition de 
certain points. Le nombre de points a apparaitre ou a disparaitre est determine par la 
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chaine de nombres poissoniens entre et L'algorithme de migration est le 
suivant : 

a) Calculer le nombre maximal iV^ de la chaine N(t) . 

b) Completer Tensemble 1 de A^^^ -N^ points et i'ensemble 2 de A^,^ -Nj 

points. 

c) Pour chaque ensemble, classer les points de 1 ^ . Tons les points 
initiaux sont classes au debut. 

d) Calculer la trajectoire de nMgration du point n de Tensemble 1 au point n 
de I'ensemble 2. (n = 1,2,.,,, ) 

e) A chaque etat t de Tensemble, supprimer les - N(t) demiers points. 
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REVENDICATIONS 

1) M6thode pour d^fonner graduellement une realisation initiale form6e 
d*aprfes des mesures ou observations et d^finissant la repartition d'un ensemble 

5 d'objets dans une zone d*un milieu h^terogene tel qu'une stmcturc g6ologique, 
generee par simulation d'un modele stochastique de type objet, les objets etant 
repartis dans la zone suivant un processus ponctuel de Poisson sous forme de points 
figuratifs avec une densite de points A(jc) qui varie en fonction de leur position (x) 
dans la zone, caract^ris^ en ce que : 

10 - on gendre une realisation d'un vecteur al^atoire uniforme selon lequel la position 
de chaque objet est definie tout en respectant la densite A(:^:) ; et 

- on modifie graduellement le vecteur aleatoire uniforme selon une procedure de 
deformation graduelle, pour obtenir la migration graduelle de chaque objet et par 
consequent le changement graduel de la repartition des objets dans la zone, 
15 jusqu*^ obtenir une realisation finale ajustee au mieux a des parametres relatifs a 

la structure du milieu tels que des parametres hydrodynamiques. 

2) Methode selon la revendication 1, caracterisSe en ce qu'on limite la 
migration des points figuratifs d*objets dans un sous-domaine de la zone en imposant 
une densite de points nulle dans la partie complementaire du sous-domaine. 

20 3) Methode selon la revendication 1 ou 2, caracterisee en ce que Ton passe 

graduellement d'une realisation contenant un prennder ensemble de AT, points k une 
realisation contenant un deuxieme ensemble de N2 points en construisant une chaine 
N{t) de nombres de Poisson entre les deux nombres iV, et N2 en utilisant la 
procedure de deformation graduelle. 
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4) M6thode selon Tune des revendications 1 3, caracterisee en ce que Ton 
modifie graduellement la taille, la forme et Torientation d'un objet au cours de sa 
migration en utilisant la procedure de deformation graduelle. 

5) Methode selon Tune des revendications 1^4, caracterisee en ce que Ton 
ajuste graduellement la densite de points en utilisant la procedure de 
deformation graduelle. 
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